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ABSTRACT 

We investigate the accuracy attainable by forthcoming space-based observations of the cosmic microwave back- 
ground (CMB) temperature and polarization anisotropy in constraining the dark energy density parameter £1q and 
equation of state wq = Pq/pq- Despite degeneracies among parameters, it is possible for high precision observa- 
tions such as those from the Wilkinson Microwave Anisotropy Probe and Planck to provide interesting information 
on the nature of the dark energy. Furthermore, we show that imposing a flat universe constraint makes it possible 
to obtain tight limits in the space of dark energy parameters even from the CMB alone. 

Subject headings: Cosmic microwave background — Cosmology: observations — Cosmology: theory — 
equation of state 



1. INTRODUCTION 

Recent cosmological observations have shed new light on the 
cosmic budget problem. There is now strong evidence from 
CMB anisotropy measurements (see e.g., Benoit et al. 2002 and 
references therein) that the universe has a total energy density 
close to critical (and therefore a flat large scale geometry) and 
that matter (either luminous or dark) can only account for about 
30% of the total density (see e.g., Turner 2002). On the other 
hand, direct evidence for cosmic acceleration from high red- 
shift type la supernovae observations (Riess et al. 1998; Perl- 
mutter et al. 1999) can be interpreted as indication for the ex- 
istence of a smooth dark energy component with equation of 
state wq = Pq/pq < -1/3, that would be responsible for the 
remaining 70% of the critical density. The nature of this dark 
energy component, however, remains a mystery. In particular, 
little is known about its equation of state. The simplest possible 
kind of dark energy is the vacuum energy (or cosmological con- 
stant), with wq = -1 independent of time. However, any plausi- 
ble scalar field from fundamental theories has a vacuum expec- 
tation value that would overdose the universe by tens of orders 
of magnitude. More general scalar fields, termed quintessence, 
may be spatially inhomogeneous, with wq ^-1 and varying in 
time. A number of strategies have been proposed to investigate 
the nature of dark energy: supernovae observations at multi- 
ple redshifts (Huterer & Turner 2001), weak lensing (Huterer 
2002), cluster counting (Haiman, Mohr & Holder 2001), red- 
shift surveys (Matsubara & Szalay 2003), Lyman-a forest (Viel 
et al. 2003), CMB anisotropy (Baccigalupi et al. 2002), as well 
as combinations of these methods (see, e.g., Kujat et al. 2002 
and references therein). 

In this work, we investigate the ability of space-based high- 
resolution CMB anisotropy observations to constrain the dark 
energy parameters. We focus on the currently underway Wilkin- 
son Microwave Anisotropy Probe (WMAP) 4 satellite mission 
and on the forthcoming Planck Surveyor 5 and we produce pro- 
jections for the model parameter errors that are attainable by 
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these experiments. We model the dark energy component as 
a minimally coupled quintessence field following tracking tra- 
jectories for an inverse power law potential (Ratra & Peebles 
1988; Wetterich 1988). Such a field is in general character- 
ized by two main features: (1) spatial inhomogeneities are on 
scales comparable to the horizon (whereas in non-minimally 
coupled theories they can generally appear on all scales [Per- 
rotta & Baccigalupi 2002]) and (2) the equation of state does 
not vary in time during the tracking regime, when quintessence 
is subdominant; in this regime the equation of state is fixed 
by the attractor solution of the Klein Gordon equation. The 
equation of state can, however, change significantly at low red- 
shifts, when quintessence is no longer subdominant, and in this 
case it tends to a cosmological constant behavior, thus having 
a present equation of state generally smaller than during the 
tracking regime. 

2. METHOD 

To quantify how well WMAP and Planck can estimate the 
dark energy parameters, we followed a Fisher information ma- 
trix approach (Fisher 1935; see also Tegmark, Taylor & Heav- 
ens 1997). The Fisher information matrix is defined as the ex- 
pectation value 

where L is the likelihood of the data and p, are the model param- 
eters. The diagonal elements of the inverse Fisher matrix give 
the minimum variance of model parameters estimated from a 
given dataset, assuming an underlying fiducial "target" model. 
This approach relies on the assumption that the likelihood is 
well approximated by a Gaussian around its peak. 

The Fisher matrix depends on the cosmological model and 
parameters adopted, as well as on the covariance matrix of the 
data, which in turn depends on the experimental setup. When 
only temperature information is considered, the Fisher matrix 
can be calculated from the theoretical power spectrum of CMB 
anisotropy C\ and its first derivatives with respect to cosmolog- 
ical parameters p,-, according to the formula: 

For an experiment sampling a fraction of sky / s k y with N fre- 
quency channels c, each with average sensitivity per pixel a c 
and a Gaussian beam response B c with FWHM angular res- 
olution 6 C , if we define w = ^ c w c = Y^S a <:9cT 2 and Bj = 
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Table 1 
Experimental Parameters 







WMAP 




Planck/LFl 




Planck/UFl 




Center frequency (GHz) 


40 


60 


90 


44 70 100 


100 


143 217 


353 


Angular resolution (FWHM, arcmin.) 


31.8 


21 


13.8 


24 14 10 


9.2 


7.1 5.0 


5.0 


a T per pixel ( x 10" 6 ) 


4.1 


9.4 


21.8 


2.7 4.7 6.6 


2.0 


2.2 4.8 


14.7 


a p per pixel (x 10" 6 ) 








3.9 6.7 9.3 




4.2 9.8 


29.8 



Note. — Sensitivities to temperature and polarization, <x and o , are relative to the average CMB temperature (2.73 K). A pixel is a square whose side is the 
FWHM extent of the beam. 



^B 2 cl w c /w, we can write (see Knox 1995): 

(AC ' ,2 *m^K, (c >^' BJ1)1 (3) 

When the polarization information is also considered (as in this 
Letter), then equation (2) has to include the full covariance ma- 
trix of the C; for all the temperature and polarization compo- 
nents: the full formalism, which we adopt in our analysis, is 
given, for example, by Zaldarriaga, Spergel & Seljak (1997). 

The WMAP and Planck experimental parameters assumed 
in our analysis are summarized in Table 1. We included in 
the analysis three frequency channels from WMAP, three from 
Planck Low Frequency Instrument (LFI) and four from Planck 
High Frequency Instrument (HFI). We fully took into account 
the Planck polarization capabilities. For both experiments, we 
conservatively assumed that foreground contamination from the 
Galactic plane leaves unobserved a symmetric strip of ±20° 
around the Galactic equator, so that the observed fraction of 
the sky is / s k y = 0.66. The mission duration required to meet 
the experimental specifications used in this paper is 2 years of 
continuous observation for WMAP and 14 months (two sky sur- 
veys) for Planck. 

We considered inflationary adiabatic cosmological models 
with eight free parameters: the present-day dark energy den- 
sity £!q and equation of state wq; the total energy density of 
the universe ft; the physical baryon and cold dark matter densi- 
ties Wb = flbh 2 and w c d m = f2 c d m /z 2 (here h is the present value of 
the Hubble parameter in units of 100 km/s/Mpc); the primordial 
spectral index of scalar perturbations n s \ the overall amplitude 
of the CMB power spectrum in units of the CQBE'-normalized 
Cio multipole; the ratio R between the tensor and scalar contri- 
bution to the CMB quadrupole. We used a single-field inflation 
consistency relations between the tensor amplitude and spectral 
index: nj = -R/6.8. We assumed no massive neutrino contri- 
bution and a negligible reionization optical depth. Note that 
£1 = Qm + ^q, where £1 M = fib + ^cdm is the total matter den- 
sity: therefore h is a dependent parameter, determined by the 
constraint: h= [(wb + w c dm)/(^ - ^Q)] 1,/2 - The choice of Q, uj\, 
and cj c dm as free parameters (rather than, for example, fib, fi c dm 
and h) has become usual practice in this sort of analysis. In 
fact, the combinations uj\, and uj c ^ m directly govern the physics 
of acoustic oscillations which defines the CMB anisotropy pat- 
tern, while il fixes the geometry of the universe and then the 
angular size of characteristic features on the CMB. For this rea- 
son, fi, Wb and w c d m are much better constrained by the CMB 
than other combinations of parameters, and are therefore a more 
suitable choice in a Fisher matrix analysis (Efstathiou & Bond 
1999). 



Table 2 
Target model 



Parameter 


Value 


wq 


-0.8 


fi Q 


0.7 


n 


1 


^b 


0.022 


Wcdm • • • • 


0.145 


n s 


1 


R 


0.1 


h 


0.746 



Note. — Cio is normalized to COBE. 



We chose our target to be the flat quintessence model which 
best fits the currently available CMB data (Baccigalupi et al. 
2002). The parameters of this model are summarized in Table 2. 
To quantify variations around the target model we numerically 
computed two-sided derivatives of the theoretical CMB angular 
power spectrum with respect to the parameters, using a step size 
which was roughly 5% of the target parameter value. All the 
theoretical spectra were computed using a modified version of 
CMBFAST (Seljak & Zaldarriaga 1996). The variation of the 
total energy density was taken into account through its effect 
on the angular diameter distance, which, for fixed Wb and u c & m , 
just results in a shift in multipole space of the CMB angular 
power spectrum. 

3. RESULTS 

The main effect of a dark energy component on the CMB 
anisotropy pattern is purely geometric. Varying the dark en- 
ergy equation of state changes the angular diameter distance (by 
changing the expansion rate of the universe), resulting in a shift 
of features in the angular power spectrum of the CMB towards 
larger angular scales (lower multipoles) as wq gets larger than 
-1. Varying the total energy density of the universe, O, also 
changes the angular diameter distance, because of the geodesic 
deviation of CMB photons from recombination to the present. 
We thus expect to observe a degeneracy between wq and O, i.e. 
a variation in the angular diameter distance due to wq can be 
compensated by an opposite variation due to £1. This is just an 
aspect of the well-known geometrical degeneracy (Bond, Efs- 
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tathiou & Tegmark 1997; Zaldarriaga, Spergel & Seljak, 1997) 
inherent in any CMB anisotropy measurement. The degen- 
eracy is not exact because different values of wq result in a 
different integrated Sachs-Wolfe contribution at large angular 
scales: however, this is precisely where the cosmic variance 
uncertainty on the CMB angular power spectrum is larger. The 
amount of degeneracy can be quantified by investigating the co- 
variance between wq and 0, obtained from the 2x2 submatrix 
of F' 1 corresponding to this pair of parameters. In Figure 1 we 
show the 68% confidence level constraints in the (f2, wq) plane 
obtained with this method. Clearly, varying wq even by a con- 
siderable amount has a much weaker effect than varying Q. As 
a result, the dark energy equation of state is poorly determined 
by CMB observations, even though the high sensitivity achiev- 
able by Planck allows one to put an upper limit to wq. On the 
other hand, the determination of Q is not very much affected by 
variations in the dark energy equation of state, because of the 
much stronger dependence of the angular diameter distance on 
£1. Note that these results do not change when CMB polariza- 
tion is included in the analysis. 

Constraints in the space of dark energy parameters (£!q, wq) 
obtained with the same technique are shown in Figure 2. Again, 
there exists a strong degeneracy between the two parameters: 
due to this, WMAP is basically unable to distinguish our target 
model from a cosmological constant case wq = -1. The situ- 
ation improves when a flat universe (£1 = 1) is assumed. This 
additional constraint partially breaks the degeneracy, reducing 
the allowed region in the dark energy parameter space. The 
improvement is dramatic for Planck: the confidence level con- 
tours get closed, enabling an accurate determination of the dark 
energy parameters. 

Our results are summarized in Table 3, where we show the 
projected lcr error bars on each parameter of our model when 
marginalizing the others. As it is well known, some cosmo- 
logical parameters are constrained to very high accuracy by 
the CMB: notably, the baryon and cold dark matter physical 
densities, which control the relative peak heights in the angu- 
lar power spectrum, can be measured to an accuracy of better 
than percent. Geometrical effects do not change this result con- 
siderably, which therefore remains valid when different dark 
energy equations of state are considered. Even in the presence 
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FIG. 1 . — The 68% confidence level constraints in the (S7, h>q) parameter 
space. The shaded regions (lighter to darker) are obtained from WMAP and 
Planck (temperature only) and from Planck (temperature and polarization). 



of a strong geometrical degeneracy, Planck is able to set in- 
teresting constraints on the dark energy equation of state: our 
target model wq = -0.8 can be distinguished at lcr from the cos- 
mological constant case wq = -1 using Planck measurements. 
Imposing the flat universe constraint 0=1 results in a decrease 
of roughly a factor of 2 in the error bars for the dark energy 
parameters wq and Oq. 

4. DISCUSSION 

In this Letter, we quantified the capability of WMAP and 
Planck to constrain dark energy cosmologies. The issue of 
the importance of CMB measurements for dark energy has re- 
cently received considerable attention and different interpreta- 
tions: for example, Baccigalupi et al. (2002) have shown that 
with a strong prior on the Hubble constant one could constrain 
the dark energy at the level of 10% even using present CMB 
data, also obtaining an indication that wq ~ -0.8. On the other 
hand, Bean & Melchiorri (2002) have shown that when that 
prior is relaxed, our knowledge of dark energy from CMB data 
alone is still poor. 

For the background cosmology, we evaluated the attainable 
precision of measurements of the dark energy abundance and 
equation of state, the amount of baryons and cold dark mat- 
ter, the Hubble constant and the cosmological curvature; for the 
early universe, we considered the perturbation normalization, 
the scalar spectral index, and the ratio between tensor and scalar 
perturbations, by only assuming single-field inflation consis- 
tency relations between the tensor amplitude and spectral index. 
We found that the performance of Planck is a factor 4 to 6 bet- 
ter than WMAP in all cases. For the tensor-to-scalar ratio, the 
inclusion of polarization for Planck allows to gain a factor 2 in 
accuracy, as expected since polarization is directly susceptible 
to tensor perturbations. 

As for the dark energy, we confirmed expected parameter de- 
generacies: indeed, the main effect of wq > -1 on the CMB is 
to decrease the distance to last scattering; however, such effect 
can be mimicked by a closed geometry (i.e. ft > 1). Inclusion 
of polarization information does not help breaking these degen- 
eracies. On the other hand, despite such drawbacks, our main 
result is that the capability of Planck to constrain the dark en- 
ergy is quite good; this statement holds in particular when we 
compare the P/ancAforecasted constraints on the dark energy 
parameters with those from dedicated experiments such as the 
Supernova/Acceleration Probe (SNAP) 6 . The latter is expected 
to reach a precision of a few percent in both dark energy abun- 
dance and equation of state. According to our results, Planck 
should perform at the level of a few percent accuracy on the 
dark energy abundance, but of only about 10% on the equation 
of state, mainly because of the degeneracy with the spatial cur- 
vature; if the latter is set to zero (i.e. if ft = 1), the equation 
of state can be determined at the level of about 5%. We also 
recall that, as it is well known, contours in the (Hq, wq) plane 
obtained from supernova measurements are roughly orthogonal 
to those obtained from the CMB, thus maximizing the accuracy 
that can be achieved by combining the two methods of obser- 
vation. Due to the precision attainable by Planck and SNAP, it 
is reasonable to expect that an accuracy of about 1 % on the de- 
termination of the dark energy parameters will be obtained by 
complementing the results of the two experiments. 

We point out that the results obtained here do not depend sig- 
nificantly on the target model assumed or on the quintessence 

6 See http://snap.lbl.gov. 
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Table 3 
Marginalized la errors 



Parameter 


WMAP 


Planck, T only 


Planck, T and P 


Awq .... 


0.54 (0.37) 


0.11 (0.050) 


0.095 (0.041) 


An Q .... 


0.22 (0.15) 


0.044 (0.020) 


0.037 (0.016) 


An 


0.0086 


0.0017 


0.0014 


Aui b .... 


0.00039 (0.00038) 


9.6 x 10" 5 (9.6 x 10" 5 ) 


8.2 x 10" 5 (8.2 x 10" 5 ) 


Aw c dm • • 


0.0041 (0.0040) 


0.0018 (0.0015) 


0.0012(0.0011) 


An s 


0.011 (0.011) 


0.0042 (0.0034) 


0.0031 (0.0027) 


AC10/C10 


0.044 (0.044) 


0.041 (0.041) 


0.020 (0.020) 


AR 


0.071 (0.071) 


0.054 (0.053) 


0.025 (0.024) 


Ah 


0.29 (0.20) 


0.061 (0.028) 


0.051 (0.022) 



Note. — Errors relative to the target model of Table 2. The number in parentheses were obtained by imposing the additional constraint Q = 1 . 



details; indeed, we repeated the same analysis by adopting an 
effective dark energy model, having a constant equation of state, 
and another target model, with equation of state wq = -1, ob- 
taining similar results. 

We conclude that the CMB is one of the most sensitive ob- 
servables to the main dark energy parameters, i.e. its abundance 
and equation of state. In fact, the recently released results from 
WMAP's first year of observation (Bennett et al. 2003) are en- 
couraging in this sense (Spergel et al. 2003): the actual parame- 
ter constraints are somewhat broader than estimated in our fore- 
cast, but they should get tighter when more data are collected. 
By combining CMB observations with those from dedicated ex- 
periments such as SNAP, we can expect to measure the equation 
of state at the level of percent. Finally we also recall that the use 
of CMB for any cosmological measurement, including dark en- 
ergy, has the great advantage that CMB perturbations are linear 
and therefore relatively simple to describe. With any other ob- 
servable from large-scale structure, we need to fully understand 
the non-linear structure formation in dark energy cosmology. 

The authors would like to thank Lloyd Knox and Paolo Na- 
toli for useful discussions. 
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